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Abstract 

The goal of this work is to introduce a local and a global interpolator in 
Jacobi-weighted spaces, with optimal order of approximation in the context of 
the p-version of finite element methods. Then, an a posteriori error indicator of 
the residual type is proposed for a model problem in two dimensions and, in the 
mathematical framework of the Jacobi-weighted spaces, the equivalence between 
the estimator and the error is obtained on appropriate weighted norm. 


1 Introduction 

In this paper we show several results concerning the two dimensional Jacobi-weighted 
spaces and we introduce a local and a global interpolator with optimal order of ap¬ 
proximation in the context of the p-version of finite element methods (FEM). Then, we 
consider a two dimensional model problem and we introduce an a posteriori error esti¬ 
mator of the residual type for the p-version of FEM, and we prove that the estimator is 
equivalent to the error on appropriate Jacobi-weighted norm up to higher order terms. 

It is well known that the development of a posteriori error indicators and adaptive 
procedures play nowadays a relevant role in the numerical solution of partial differen¬ 
tial equations. In contrast to the case of h refinement, it seems to be an open question 
whether uniform reliability and efficiency can be achieved for an hp a posteriori esti¬ 
mator of the residual type even in simple problems. 

In the one dimensional case, analysis for a posteriori error indicators based on the 
residuals for the p and hp of FEM is well known [ Sch98l GB861IDH07] . More precisely, 
in 1DH07 ] the authors obtain an error estimator of the residual type for the Poisson 
Problem and prove that the H 1 norm of the error is equivalent to the error estimator 
up to higher order terms. Moreover, they propose an adaptive algorithm and, since the 
error estimator is reliable and efficient, they prove that the algorithm leads to a uniform 
monotone decrease of the energy error in every step. It is important to point out that 
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these kind of results have not been established in high dimensions, and the techniques 
used for one-dimensional analysis can not be applied to higher dimensions. 

In the two dimensional case, to the best of the authors’ knowledge, the error estima¬ 
tors of the residual type present in the literature for the p and hp of FEM are equivalent 
to the error with constants depending on p (see [APRS111 APRS12[ 1APS14 . MWOlj 
and the references therein). In particular, in |MW01] the authors obtain an error es¬ 
timator of the residual type for the Poisson Problem in two dimensions and, using 
optimal weighted inverse estimations, prove that the error is equivalent to the indicator 
and propose an hp strategy based on a predictor of the error in each element of the 
mesh. Using this error indicators, a particular algorithm is proposed in . 111)1 11 and its 
convergence is reached assuming a data saturation which, due to the constant p de¬ 
pendence, becomes more restrictive for increasing polynomial degrees. In [ AS97 , AS98| 
the authors proposed an hp refinement strategy in which in every step and for every 
element they decided to do h adaptivity or p adaptivity based in the local regularity of 


the solutions. On the other hand, in 


p-robust equilibrated residual error 


estimates are obtained for the Poisson problem and Elasticity problem respectively. 

In recent decades, the Jacobi-weighted Sobolev spaces have received increasing at¬ 
tention for the approximation theory of the p (and hp) version of the finite element 
methods. These spaces seem to be the appropriate functional spaces for a priori error 
analysis and play a crucial role in the analysis of the a posteriori error estimations 
(see |AP021 IBG001 IBG02al IBG02bl IGB13] , and the references therein). Indeed, the a 
priori error analysis and optimal convergence for the p-version of FEM in this context 
have been studied by several authors (see, for instance, |BG02al [BG02bl FGS07j ). More 
recently, increasing attention to this framework has developed because of the need for 
optimal a posteriori error estimates [ BG101 lGuo05j . Motivated by the results obtained 
by [DII07 in the one dimensional case and the a posteriori error analysis given by 
[ Guo05j . in this paper we analyze the a priori and a posteriori approximation theory for 
the p-version in the mathematical framework of the Jacobi-weighted Sobolev spaces. 
In fact, we present several results concerning the interpolation theory for functions in 
Jacobi-weighted Sobolev spaces and, for the two dimensional Poison model problem, we 
develop an a posteriori error estimator of the residual type for the p-version of FEM. 
We analyze the equivalence of this estimator with the error in a Jacobi-weighted norm 
and we prove quasi-optimal global reliability and local efficiency estimates, both up to 
higher order terms. As far as we know, simultaneous reliability and efficiency estimates, 
both with constants independent of the polynomial degree, have not been proved yet for 
any a posteriori error estimator for p finite element methods in the two dimensional case 
and our estimates may be (see [ Guo05l IGB13] ) the best result that one can expected 
for error estimator based on residual in two dimensions. 

The rest of the paper is organized as follows. In Section [2] we show several results 
concerning the Jacobi-weighted spaces. In Section [3] we present the p-approximation 
theory and the interpolations error. In Section [I] we consider a two dimensional model 
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problem and we introduce the a posteriori error estimator and we prove its equivalence 
with the error in a Jacobi-weighted norm. 


2 Jacobi-weighted Sobolev spaces 


Let Q = (—1, l) 2 the reference domain in R 2 . For % — 1,2 let be > —1, cq > 0 
integer, (3 = (/Jl,/^) and a = ( cti,a. 2 ). We define the weighted function Wp tQ in Q as 
follows 

W Pta (x,y) = (l- X 2 )P 1+ai (l-y 2 ) p2+a2 . 

If a = 0 we note Wp = Wp$. 

For a function u G C°°(Q) and k > 0 integer, we define the following norm: 

= ^ U \ W/3, a - 

|a|</c Q 

The weighted Sobolev space H k ^{Q) is defined as the closure of the C°°(Q ) functions 
with this norm (see, for example, |BG02aj ). i.e., 

H k ' p {Q) = . 

With \u\ H k t piQ) we denote the seminorms 

\u\hW{q) = £ / \d a u\ 2 Wp, a . 

\a\=k^ Q 

Let be an open polygonal domain in R 2 , T an admissible partition of Q in paral¬ 
lelograms. For any K G T, let F : Q —> K be an affine transformation and u G C°°(K), 
then u = u o F G We define 

— ll^ll 5 (1) 

and 

II 112 \ ^ || 112 

\\ U \\H k ’P(T) = INI H k ’P(K)' 

KeT 

We observe that the norm depends on the partition that we are considering. Then, the 
Jacobi-weighted spaces for T, a partition of 0, is defined as 

H k 'P(T) = C'oo(Q) IMI ^( r) 

From now on, we consider the case /3j = f3 > — 1 for i = 1,2. 
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In what follows we present the Jacobi projection and we enunciate its properties 
(see [GS07] and the references therein). In order to do that we need to introduce the 
Jacobi polynomials in one dimension (for details see [ Guo09] ). 

Let / = (—1,1), (3 > — 1 and let p be a polynomial degree. For iGl, let (x) be 
the Jacobi polynomial of degree p, i.e., 

0 (1 - x 2 )~P d p { 1 - x 2 f +p 

p X ~ 2Pp\ dx p 


It is well known that the Jacobi polynomials Jp(x ) are orthogonal with the Jacobi 
weight Wfj(x) = (1 — x 2 ) 0 , i. e., 


J 4(x)X m (x)W f (x) 


7p, p = m 

0, p m 


with 

_ 3 2 2 h+ 1 T 2 (p + (3 + l) 

lp ~ ( 2 p + 2(3 + l)r(p + l)r(p + 2/3 + 1) ’ 

where T denotes the well known function Gamma given by T(z) = J 0 °° t z_1 e~ l dt. 
Let Jp k {x) = -£jrjP(x), then for 0 < k < p we have 



-k ^(P + 2/3 + k + 1) 0+k 
T(p + 2/3+1) p ~ k[ ) 


which are orthogonal with the Jacobi weight Wp + k(x ); 


J rA x ) J ™,k( x ) W P+k (x) 


7 p,ki p = m> k 
0 , otherwise 


( 2 ) 


with 

_ 2 2/3+1 r(p + 2(3+ k+ l)T 2 (p + P + 1) 
lp ' k ~ (2 p + 2(3 + l)r(p + 1 - k)T 2 (p + 2(3 + 1)' 

We note that if k = 0 we obtain 7 = 7 ^. 

Now, we will enunciate two important properties that provide us an estimation for 
the constant 7 ^ k and 7 ^. For this purpose we need some previous lemmas. 

The following lemma (which can be found, for example, in page 427 of |Lan99j l 
gives a well known estimation for the function F. 


Lemma 2.1. For real x and x —> +00 the following applies 

r(x) - x x ~ 1/2 e~ x VZn 


where ~ means the quotient of the left side by the right side tents to 1 as x —» +00. 
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Thus, we have the following useful results. 

Corollary 2.1. a) For any a Gl, r(n+a) ~ (n+a) n+a 1 / 2 e ( n+Q ^v / 2vr, for n —» 
+oo 

b) For any aeR, lim „^ +00 ^^*2 = 1, 

c) Given ap > — 1, for any a such that — 1 < a < a 0 , there exists positive constants 
A and B (depending on aio but independent of a) such that 

A < T j\ + ^ < B Vn G N 
~ Y(n)n a ~ 

Therefore, the following estimate holds. 


Lemma 2.2. Let —l</3</3o,k>0 integer and p a polynomial degree. There exist 
positive constants A = A((3 0 ) and B = B(/3 q), independent of p and (3, such that 

Aiy l < -rh, < Bp 2k Vp > k. 

Proof. If k = 0 we take A — B — 1. Suppose k > 1. 

7 p, k _ 2 2 P +1 r(p + 2(3 + k + l)r\p +(3+1) (2p + 2(3 + l)T(p+l)T(p + 2(3 + 1) 1 

^ p 2k ~ (2p + 2(3 + l)T(p+l-k)T 2 (p + 2(3 + 1) 2 2 h+ 1 T 2 (p + (3 + 1) p 2k 

T(p + 2(3 + k + l)T(p + 1) 

T(p + 1 — k)T(p + 2(3 + 1 )p 2k 

T((p + l-k) + (2P + 2k)) T {p)p 2 P +1 rp + l-kyh+2k 
~ T(p + 1 - k){p + 1 - k) 2/3+2k T(p + 2/3 + 1) V p ) 

= (/)(//)(///). 

In (/) and (II) we apply Corollary 12.11 c). for (III) we observe that 


p + l-k 
inn - = 1 , 


p—kOO 


P 


and, since p+1 p k > 0 for all p > k, we get 

p + l-k 

ci < --< c 2 , Vp > k, 

p 

for positive constants C\ and c 2 and the proof concludes. 


□ 


Lemma 2.3. For —l<(3</3o and p > 1 a polynomial degree there exist positive 
constants A = A((3 0 ) and B = B((3 0 ), independent of p and (3, such that 

Ap~ x < 7p < Bp 1 . 
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Proof. 


_ 3 2 2 P +1 r( P + /3 +1) 2 

lp “ ( 2 p + 2 p + i)r(p + i)r(p + 2/3 + 1) 

2 2 p+i , r (p + p + i)y p 2 P +l v(p) 

~ ( 2 p + 2/3 + 1) l Y(p)p^ 1 ) r(p + 2/3 + 1 ) 


By Corollary 12.11 c) there exist positive constants c\ and c 2 , independent of p and /3, 
such that 


/ T (p + f3 + 1) 

Cl ~ v r ^)^ 1 


2 p 2 ^ +l r(p) 

r(p + 2/3 + 1 ) 


< c 2 ,Vp > 1,V - 1 < /3 < /3 0 


and the result holds. □ 

For any /5 > —1, k > 0 and u G H k ’P(Q), we have the Jacobi-Fourier expansion 
(see, for instance, |GS07j ) 


with 


u(x,y) = Y c i,j J i( x ) J j(y) 

i,j =0 


u(x)J?(x)j!l(y)Wp(x,y). 


(3) 


C i,j 8 8 . 
lil j JQ 


Using the orthogonality of the Jacobi polynomials ([2]) we have that 

OO p 

Y iGjf 'ihj and / \d a u\ 2 W 0ta = Y ( 4 ) 


a \\H°’P(Q) ~ 

i,j =0 

Then, for —1 < (3 < j3o, by Lemma [+2l we deduce that 
12 


l>OL\ ,j>CX 2 


m 


, „ -V V \r- -iV* V 9 ~ V V | r . .|2^+2 qi -2a 2 

H k ’P(Q) 2-^ 2-^ +*+ 7i,ai^j,a 2 — 2^ 2^ +*+ 3 ’ 

|a|=fc i>ai J>a 2 |a|=fc i>ai J>a 2 


where A ~ B means ciU < A < c 2 i? with positive constants c\ and c 2 independent of 

/+ 

For p > 0 we define Q P (Q) the set of all polynomials of degree less or equal than p 
in each variable in Q. The Jacobi projection of u in Q P (Q) is 

p 

n %u(x, y) = Y Gj J f (■ x )Jj(y )• (5) 

i,j=0 
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Lemma 2.4. Let — 1 < (3 < (3o, then 


|^(-l)| = |^(l)|~ F ^i rT y( P +l/ 

with constants depending on (3 q but independent of p and (3. 

Proof. By equation (2.2) of |! GW04| we know that 

iy<-i)! = iy<i)i. aud j|( i ) = r ^t)f ( ^+ 1) - (6) 

Then, we can write 

/an i = r( P + /j + i) (p + iy 
’ r(p + i)( P + i)ar(/3 +1)’ 

and the result follows from Corollary 12 .II cl □ 

The following theorem gives approximation properties of the Jacobi projection. 

Theorem 2.1. Let I = (—1,1), Q the reference domain in M 2 , —1 < (3 < (3o, and 
u G H 1,/3 (Q). Let p > 1 a polynomial degree, II G Q P (Q) as in Then there exist 
a positive constant C = C((3<f), independent of u, [3 and p, such that 

||n — IIpit||^o,/ 3 (Q) < C(p + 1) 1 \u\ H i,p(tQy (7) 

If, in addition, u G C°(Q) and j3 < — 1/2 then 

\\{u-U^u)(±l,y)\\ H0 ,i> {I) < ^ (p+ 1)~ 1/2 \u\ h i,p { q) , 

C 

\\(u-U^u)(x,±l)\\ H o,p {I) < + (p+ iy l/2 \u\ H i,p {Q) , 

and if, in addition, (3 < —1/2 then 

|(u-Il£u)(y)| < _ 2 ^ r — (p + lf +l/2 \u\ H i, P{Q) , V V vertex of Q. (9) 

Proof. The proof of (JTj) is given in [ GS07] with a constant that could depend on (3. 
However, following the steps of that proof, we observe that a positive constant C can 
be chosen independent of (3. Now, we prove ((HJ) , the bound in the edges of Q for 
smooth functions. We carry out the case ||(u — Ii^u)(x, —l)||#o,/?(/), the other cases can 
be obtained analogously. Since u G C°(Q) we can write 

(i/-nju)(x,-l) = ( ) c hi J i( x ) J j (-!) 

2 >p+lj>p+l i>p+l,j<p+l i<p+l,j>p-\-l 

= v dGw+ X if J/W+E 

z>p +1 2 >p+l i<p +1 
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where 


bf ] = Y djJji- 1 ), i > p + 1 


j>p +1 


b ? ] = c h3 J j(~ l )i i>P + 1 


j<p +1 


h i - ^2 c ^ J f( _1 )> *<p+i- 

j>p+i 

It is well known that, if / : \p, +oo] —)■ M is a non negative, decreasing and integrable 
function then, 

00 />oo 

Y /( n ) - / f ( x ) dx ■ ( l0 ) 

i=p +1 J V 

Therefore, by Lemma [2.41 Holder inequality, Lemma [2.31 and (HUD we get, 

id‘'i 2 < ( E MIJ?(-i)l) 2 

j>p +1 

r(/j + 1)2 ( I] |cMl(j + l)^7f) 1/2 (7f)" 1/2 JJ _1 ) 2 


< 


i>p+i 


< 


C 


L(^ + l ) 2 


j>p +1 


< 


(7 


rCd + i ) 2 


(E i^i 2 7^- 2 )( E^ + i) 2 "^)- 1 ^- 2 ) 

3>P +1 

(E 

i>p+i 


fclTEK E r I+w ) 


< 


<7 


j>p +1 
oo 


< 


r(/? + i) 5 

(7 

r^Tij’ 2 


_ /*oo 

(EM 2 77 (/ p - 1+2 "<7 


i>p+i 


p 2 ' E Mhb ' 2 


j>p+i 


Then, 


^IV < r(/ ^ 1)2 p 2/j I] hjl 2 7fiV 


z>p+l 


i>p+l 




< 




n 2/3 


r(/? + 1 ) 


2 P~" J2 

i>0,j>l 


C 


12 8 -2 8 ^ ^ 2/31 | 2 

,il IjJ 7 i < Y(f) + Vf p \ u \h^{qy 











Analogously, 

if>! 21 r < (v | Ci jjf ( -i)|) 2 

j<p +1 


Thus, 


< 


< 


r(/ 3 ^ 1 ) 2 (p +1 ) 2 ( E l c M'IO' + 1 ) / 3 (^) 1 / 2 ( 7 f) 1/2 *)‘ 


< 


E IcylVhK E W + iW)- 1 ) 

' ' j<p +1 j-cp+l 

^ ryrnF (!,+ir2 ( £ wV'-H E u + u^ 1 ) 


i<p+i 


i<P+i 


< 


r^ny(f' + 1 )' 2 (E i l^l 2 -'A)(p + 1 ) 


E "i 2 h 


J 2 1 7^r x ')ii 2 


bpA-rf 


2 >p+l 


< 


Similarly, 


As a consequence, 


\h°’P(I) — E I fc ! ; 

2 >p+l 

r( «+ n. ^‘ 1' 1 E MV’V 

^ ' i>p+l,j<p+l 

- r(fl+ E l^l'AV 

c 

- r( / 0+i) 2 ^ +1 )~ 1 l w l^w)‘ 


E ^ ] '^f 1 ( X )ll/ 7 °^(/) ^ _|_ ]_) 2 -P VltfCAQ)- 


||(m - nJu)|| H o,p (7) < ^ (p + 1) 1/2 \u\ H i,P iQ) 

and © holds. 

Finally, we prove (|9|) for V — (-1,-1), the same arguments can be used for the 
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other vertices of Q. 

|(u-nj«)(-i,-i)| = ( E + E 


+ E 

i>p+l,j>p+l i>p+l,j<p+l i<p+l }< 7>p+l 

<( E + E + E )i<yi^(-i)HJ?(-i)i 

i>p+l }< j>p+l i>p+l,j<p+l i<p+l,j>p+l 

r(/3 + !) 2 ( XI + XI + X )|cij|(i+ 1 ) /3 (j + 1 ) /3 


< 


2>p+lj>p+l i>p+l,j<p+l £<p+l,j>pH-l 


c 


;(I + II + III). 


F(/3+ l ) 2 

At first we compute II and III. Due to Holder inequality and Lemma [2.31 we obtain 

a <( E m 2 7?A 2 ) 1/2 ( E (i+i^w+inTfr 1 ^)- 1 *- 2 ) 172 

*>p+hi<p+i *>p+hi<p+i 

<c( e mY 7?; 2 ) 1/2 ( E (<+i) w o+i) w <- 2 ) ,/2 


*>p+ij<p+i 


*>P+1J<P+1 


1/2 


<c( e i^iVifi 2 ) 1/2 ( E (<+i) w <- 2 ) ,/2 

i>p+lj<p+l *>P+1J<P+1 

<C( E l^l 2 T?Tfi 2 ) 1/2 ((P+l) Em 2 '’ -1 ) 172 

i>p+l,j<p+l i>p+l 

Now, from (flUl) we follow that 

II<C( Y \ c i,j\ 2 lilj i2 ) 12 {(p+ 1 ) / x _1+2/3 dx) 

*>p+ij<p+i p 

<c( e i^iV7h 2 ) 1/2 (( P +i)P 2,, ) 1/2 

*>p+ij<p+i 

<c( E l^lV7h 2 ) 1/2 (P 2,i+1 ) I/2 <C/ +1/2 |“l»«W l - 

*>p+lj<p+l 

For /// we proceed analogously but changing the roles of i and j, and we can 
conclude that 

III<CpP +1 ' 2 \u\ H i, HQ) . 

In what follows we compute I. 

/<( y idpIVt f*i) 1/2 ( X ( i + i ) 20 u+ i ) 2P (iiy 1 ( r rj)~ lrl r i ) 1/2 

2>p+lj>p+l 

<c( E i^i 2 7?7fv) 1/2 ( E^ E «+ 1 ) 2 ^) 1/2 


£>p+l,j’>p+l 


i>p+l j>p+l 


<C( X] I c m|V7/D) 7 ( / X 2 ?dx y W dy) 7 . 

i>p+lj>p+l ^ ^ 
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Since (3 < — 1/2 the integrals converge, and it follows that 


/ < 

c 

( E M 2 

-1-2/3 



< 

C 

p w ( E 

1 

1 

to 


£>p+l,j>pH-l 

< 

c 

P W+ 1 \u\m.P( Q ). 

-1 - 2/3 




d/2 


( V |c i ,|V 7 fi 2 + V 

z>p+l,jf>p+l 


Therefore, 

|(«-nj«)(-i,-i)| < ( _ 1 _ 2/3^r(/3 + i) 2 ^ + 1 / 2 I m I 

and the proof concludes. □ 

Let ih be a parallelogram of M 2 and let Fk : Q K be an affine transformation. 
For p > 0 we define 


Q p (/i) = {M|no^e Q P (Q)}. 

Now, from (jT]) and Theorem 12.11 we have the following result. 

Corollary 2.2. Let K be a parallelogram of Ml 2 , —1 < f3 < f3 0 , and u G H 1,/3 (K). 
Let p be a polynomial degree, there exist Li^ K u G Q P (K ) and a positive constant C 
independent of p, j3, and u such that 

ll M — ^p,k u \\h°’P(k) < C{p + 1 ) l \u\ H i,fi( K p ( 11 ) 

if, in addition, u G C°(K ) and f3 < —1/2 £/mn 

C 

W u ~ U p,K u \\HM(-y) < + ^ (P + 1 )~ 1,2 \u\h^P(k) V 7 edge of K, (12) 

if, in addition, [3 < —1/2 then 

l(« - U p u )( v ) I < ^_ 1 _ 2/3)T{3 + l ) 2 (P + 1 ) 1 W Mh 1 -'W V ^ vertea; °f K - ( 13 ) 

3 ^"Interpolation of smooth functions 

The goal of this section is to introduce interpolation operators which are suitable to 
obtain a residual error estimation in the mathematical framework of Jacobi-weighted 
Sobolev spaces of the p -version of Finite Element Methods. 
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Let fl be an open polygonal domain in M 2 , T an admissible partition of Q in paral¬ 
lelograms. Let p be a polynomial degree, we denote 

S P (T) = {ue C°(n) | u\ K g Q P (K ) VA g T}, 

Sq(T) = {ue C' 0 °(O) I ti|* g Q P (A) VA' g T}. 

We choose a polynomial degree px for each K G T and we note p = (px) the 
vector of polynomials degrees. We assume that the polynomials degrees of neighboring 
elements are comparable, i.e., there exists a positive constant C such that 

Pk < Cp K ' K, K' G T with A fl K' ^ 0. 

We introduce the following notation 

S P {T) = {u € C°(Si) | u\ K £ Q pk (K)}, 

SS(T) = {u £ C 0 °(!J) | u\ K £ Q pk (K)}, 1 

and 

£ = {all edges e in T}, 

£° = sn n°, 

Af = {all vertices V in T}. 

For V G A/” we denote 

CoV = (J{A|A G T and Afli'/ 0}, 

7y = T|o^, 

= min{p x |fo G A"}, 

£y = {all edges e of 5 such that V is an endpoint of e}. 

For any K G T or e G £ we define 

= (J{A'| K' G T and K' n A ^ 0}, 
w e = (J{A|A G T and KHe^ 0}, 
p e = mm.{px\e edge of A}. 

Let K G T and let Fx : Q —> A' be an affine transformation, for a function ?/ in A 
we denote u = u o A^. Let e = I x {—1}), then for any e G £ let A e : e —>■ e be an affine 
transformation then, for a function u in e we denote u = u o F e 

In order to introduce the local and the global interpolation operator we need some 
previous lemmas. 
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Lemma 3.1. Let I = (—1,1), —l</3</3o and p a polynomial degree, there exists 
g £ V P (I ) such that g{ 1) = 0, g(—1 ) = 1 and 

||^||i?°^(/) < CT([3 + l)(p + 1 ) (1+/3 \ 


where the constant C = C(j3o ) is independent of p and (3. 

Proof. Let be the Jacobi polynomial of degree p — 1 which is orthogonal with 

the weight (1 — x) /3+2 (l + x)@ (see [Guo09j for details), we define 

(1 

aw 2 jyf'»(-i) ’ 

then g G V P (I), g{- 1 ) = 1 , g( 1 ) = 0 and 


2 


y+fA-i)) 2 J‘ 

4(4d 2 '«(-l )) 2 l 

1 


(1 - x) 2 (J^ 2,/3) (x)) 2 (1 - x^l + xY 
(Jp^ 2 ,/ 3 ) (x)) 2 (l - x) 2+/3 (l + x)^ 


J + 2,/3 

2 'p—1 


As it has been shown in [Guo09] equation ( 2 . 2 ) 


4-T 2 ,/ 3 ) (- 1 ) = 


(-i)p~ 1 r(p +/?) 
(p-ijjr^ + i)’ 


and from Lemma [2.41 we have that 

dTA-i) ^ 


pr 


r(/? + i )■ 


Then, 

l^llk+r) < CT(/3 + 1)V 2/3 7^ 1 2,/3 , 

and from equation (2.7) of |Guo09] we know that 

e+w = 2 w r(p + /3 + i)r(p + /3) 

7p_1 2 p + 2/1 + 1 r(p)T(p + 2/3 + 2 ) 

2 2/3 + 3 r(p +/3 + 1 ) r(p +/?) r(p)p 2 ^+ 2 1 

2p +2/9 + 1 T(p)p^ +1 T(p)pd T(p + 2f3 + 2)p 
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Now, by Corollary 12.11 c) 


r(p + /?+1) 

i»p' 5+1 _ ’ 


r(p + ffl 

r( P )pP - ’ 


and 


t(p)p W2 r 
r(p + 2 P + 2 )~ ’ 


and hence, 

/ 3 + 2,/3 , ^ _2 

7p-i < Cp , 

and we conclude that 


Il9llt».<>(/) < CTOS + 1) V^jT 2 = CT(/? + ifp-V+D 


as claimed. 


□ 


Corollary 3.1. Let K be a parallelogram and Vi, V 2] V 2 and V 4 its vertices. Let —1 < 
j3 < /3 q and p be a polynomial degree. There exist a function and a positive constant 
C = C(j3o) independent of p and ft such that 

l) UAVj) — $ij (takes the value 1 in Vi and 0 in the others vertices), 

n) (k,i e Q P (K) , 

in) Mk,i\\h°, 0 (k) < CT(/3 + 1 ) 2 (p + l)" 2 - 2 / 3 and 

iv) ||^,dl^(e) < CT(/3 + l)(p + V e edge of K. 

Proof. Let Q = (—1, l) 2 the reference rectangle, p a polynomial degree and g as in the 
Lemma [3711 then, the function G(x,y ) = g(x)g(y) is in Q P (Q) and satisfies 

G(— 1, — 1) = 1 and takes the value 0 in the others vertices of Q, 

I|G|| h ».»( Q) < CT09 + l) 2 (p + l)- 2(1+ », 
l|G|| W(e) <CT(/3 + l)(p+l)-< 1+ » Ve edge of Q. 

Let Fxy : Q —> K be an affine transformation such that —1) = Vi then, 

(,k,i = Go F K y t satisfies i)-iv) and the proof concludes. □ 

Corollary 3.2. Let K be a parallelogram and e an edge of K. Let —1 < ft < (3q and 
p be a polynomial degree and w G V p (e) such that w\g e = 0 (i.e. w(V) = 0 for all V 
vertex of e). There exist if an extension of w to K and a positive constant C = C(j3o), 
independent of p and ft, such that 

1 ) if £ Qp(K) 

a) if \ e = w, 

Hi) if\dK\e — 0 and 
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iv) 11 "011 < C'r(/3 + l)(p + 1) ^ 1+/5 ^||M;||^o,/3( e ). 

Proof. Given Q = (—1, l) 2 , the reference rectangle and e = (—1,1) x { — 1} an edge of 
Q, we consider F Ke : Q —>• K an affine transformation such that F Ke {e) = e. Hence, 
we define w G V p (e) as w = w o e and 

0002/) = ^00 ~ l )g{y) 

where g is the polynomial of degree p introduced in Le mm a [3.11 then, G Q P (Q) and 
satisfies 


0|e = W, 

01 dQ\e = 0 , 

ll0ll-ff°^(Q) — CT(/3 + l)(p + 1) ^ 1+/3 - ) ||'u)||/fo,^(g). 

Hence, the proof concludes by defining 

0 = 0 ° Fx]e 


□ 

For each V G A/", in the following theorem we introduce a local operator Iy : 
//q ,/ 3 (H) D C°(ojv) —> S p (7y) and we present some local error estimates. 

Theorem 3.1. Let —1 < 0 < —1/2 and p a polynomial degree. For each V E Af and 
u G Hq’P(LI) fl C°(ojy) there exist an operator Iy : n C°(ujv) —> S p (7y) and a 

positive constant C independent of p, 0 and u such that 

\\u — Iyu\\ H o,p(K) < _i — 2 /3 ^P ^ ( 3 ^ 2 +/ 3 ) I m I^ 1 ’ /j (Tv) e 7v, 

C 

\\u — Iyu\\ H o,p( e ) < + 1 ) 1 ' / 2 |' u l^ 1 -^(rv) e £y- 

7/7 G (5 \ £°) is srtc/i 7 C <9ov t/ien /yu | 7 = m | 7 = 0. 

Proof. For each K G 7y we consider K u as in Theorem 12.21 Now, we will define 
f>K such that 4>k(V) — U {V) for ah V vertex of K. Indeed, given If, V 2 , V 3 and V 4 the 
vertexes of K, we define the polynomial fx of degree p as following: 

4 

&K = U p,K U + J2( U ~ U p,K U )( V l)£K,l 
1=1 

with fx.i the function defined in Corollary 13.11 
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Therefore, using (fTTj) , (jT3]l and Corollary 13.11 iii) we find that 


\ u — ^kWhO’P^K) < \\ u — H-p,K U \\h 0 >P(K) + £l(“- n pV)M)lll&dl H°’P(K) 


1=1 


<C(p+l) \u\ H i,p^ K ) + C 


< 


c 


ip + 1) 


-( 3 / 2 + 13 ) 


(-1-2/3) 

u \h 1 >P(K)i 


(p+1) 


— (S/2+/3) 


U\h^P{K) 


(—1 — 2/3) 

similarly, for any e edge of K , from (jT2|) . (TT3j) and Corollary 13.11 iv). we have that 

C 


\ u - 0rf||^o./3(e) < 


(—1 — 2/5)T(/3 + 1) 


(p+ 1 ) l ^ 2 \u\ H i,f)( K y 


Let e E S° such that e C 7y D (cay) 0 and let Ad and K 2 be the elements of Tv that 
share the edge e, we consider 


w = ((p Kl - 4>k 2 ) \ e e V P (e). 


Observe that 


H|n-0,/3( e ) < ||w - + \\ u - 0i^ 2 ll^°^(e) 

c 

(p + 1 ) 1 ^ 2 (l M lir 1 .5(_ft' 1 ) + \u\h^P(k 2 ))- 


< 


■1 — 2/3)T(/3 + 1 ) 

Hence, let if; G <2 p (A'i) an extension of w to K\ as in Corollary 13.21 then. 

i>\e = W, 

fp\dK!\e = 0 , 

||' ! / , ||h°./ 3 (A'i) < CT(/3 + 1 )(p + 1) (1+/ ^ ||'ia||#o, / 3 ( e ). 

We define 4>k 1 = </>Ai — V' and we n °f e that this function satisfies: 

0Ai|e = 4>K 2 \e, 

^KildK^e = 0Ai|aAi\e, 

and therefore we have the following error estimates 
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< 


c 


< 


(_1 _ 2/3)r(/3 + 1) 

c 


(p+1) ^ 3 / 2+ ^|tt| i? i, / 3( K i) + CT(/3 + l)(p + 1) (1+/3) |Mltf°^(e) 


(p + 1) 


(-1 - 2/3)r(/3 + 1) 

+ CT(/3 + l)Go + l)- (1+ ^ 

c 


-( 3 / 2 +/ 3 ) 


c 


U \H 1 ’P(K 1 ) 


< 


-1 - 2/3 


( — 1 — 2/3)r(/3 + 1 ) 

(P + 1) (|M | h 1 ’! 3 (Ki) + 


(. P + 1) 1 ^ 2 (I' u I// 1 ./ 3 (Ki) + \ u \h 1 ’P(K 2 )) 


and 


\\ u ~ 0Ki||ff°.0(e) < ||« — <j)K! HifO./3(e) + ||V’l|if°>0(e) 

= ||w - 4>Ki§HW(e) + \\ w \\H°’P(e) 

~ (-1 - 2 /3)T(/3 + 1) ( P + ^ ' ^ Rl ’^ K i) 

C 

+ (_i _ 20 )T{P + 1) ^ + 1 ) _1/2 (I u I h1 ^(^i) + \ u \h^(k 2 )) 

C 

- (-1 - 2/3)r(/3 -h 1)^ + 1 ) _1/2 (l u l ffl ’ /5 ( Jf i) + 

Then, we can modify 0^1 to 0^, and repeat this process in each e E 8 ° such that 
e C Ty n (u)v)°■ 


If 7 E (8 \ 8 °) is such 7 C dcoy, we consider K E Tv such that 7 E <9ib. Let 
w = we observe that w — u — 0 in < 97 . Now, let 0 be as in Corollary 13.21 an 

extension of w to K and 4>k = 4>k — 0 then, 

4>k | 7 = 0, 

(pK | dK\y = 0iC0K\7, 


and 


l M - 0rdl^°>/ 3 (Ap < ll M - (PkWhO’P/k) + W^Who^k) 

-(p+1) (3 / 2+/5 )|m| h i,/ 3( K ) + CT(0 + l)(p + 1) 


< 


< 


(- 1 - 20)100 + 1 ) 
C 


(- 1 - 20 )T (0 + 1 ) 


(p+ 1) 


-( 3 / 2 + 0 ) 


W|_f/ 1 ./3(K) 


+ CT(0 + l)(p+ 1) 


-(!+/ 3 ) / 1 |, ; _ 


u ~ 4 >k\\h°’P(ci) + ll^ll j 













since ||w||#o,/ 3 ( 7 ) = 0 because u | 7 = 0 we have that 

C 


\ u - 4 >k\\h°’P(K) < 


(p + 1 ) 


— ( 3 / 2 +/ 3 ) 


< 


(_1 _ 2/3)r(/3 + 1) 

CT(^ + l)(p+l)-M||«-^||^ (7) 

-£-(p+1 )-(3/2+«- ' 


u I if 1 ./’(if) 


+ 


1 - 2/3)T(/3 + 1) ^ 1 ' UlHl ^ K) 

C (p+ l)” (1+/3) (p + l)~ 1/2 \u\ H i'P(K) 


< 


-1 - 2/3) 

—-- (p+ l)-( 3 / 2 +« 


—1 — 2/3 




We can modify to <f>K and repeat the process for all 7 such that 7 G (£ \ £°) and 
7 C <9cjy. Thus, the operator Iyu\x = satisfies all the requirements. □ 

Now, we are in conditions to introduce a global operator Iu G (T) which satisfies 
the following error estimates, 

Theorem 3.2. Let —1 < /3 < —1/2 and u G D C°(0) then, there exist Iu G 

Sq (T) and a positive constant C such that 


\u ~ Iu\\ H o,fp K ) < 


C 


\ u ~ I u \\H°>P(e) < 


-1-2/3 
C 


-( 3 / 2 + 0 ) 

Pk 


-1-2/3 

T/ie constant C is independent of u, /3 and p 


u \h 1 >p{t\u k ) V K G T, 
■P e 1 / 2 l M lni./3(ru e ) V e G £. 


(15) 

(16) 


Proof. A fundamental property of the space 5 q(T) (see, for example, |Mel05| ) is that 
we can identify “nodal shape functions” that form a partition of unity, i.e., for each 
vertex V G Af, we can find a function tp v £ S 1 ( T) such that 

</v|n\ wv = 0 , ^ <py = 1 and sup \(pv(x)\ < 1 . 


veAf 




We consider IyU G 5 ,p ' / 1 (7y) as in Theorem 13.II with p = py — 1 and we dehne 

Iu= ^ PvIyU. 


vga r 


It is clear that Iu G Cq(Q) and 


vgk 
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since py <Pk, for all K such that V E K, then Iu\k E Q Pk (K) and 


u ~ I u \\H°’P(K) 


y pvu - 'y <p v iv u \\ H °’P(K) 

VeAf VeAf 


'y, Tv{u — IyU)\\ H o,fp K ) < \\ u — Iv U \\h°’P(K) 

VeA( VeK 


< 


c 

-1-2/3 


E -( 3 / 2 + 0 ), | 

Pv \ u \m>p (t v )i 

VeK 


since polynomial degrees of neighboring elements are comparable then p v ^ 2+ ^ < 
Cp~/ /2+ P\ therefore 


II T II ^ ° —(3/2+0) | | 

If _ Iu\\hw(k) < _ 1 _ 2 pPK \ u \h^{t\u, k p 

and the hrst estimate holds. Now, let e G £ we have that 

ll M - -MltfQFe) = II y^ V yU - y, / 9 V'^V'' U llrr°’ /3 (e) < ~y \\ u - Iv u \\ n°./3(e) 

V£Ar VGAr vee 

^ C ST' —r/21 i 

- _ 23 y pv l M l 

^ Vee 

since polynomial degrees of neighboring elements are comparable then pf 1 ^ 2 < Cpf 1 ^ 2 , 
therefore 

C 

ll M — -^ M lln°./ 3 (e) < —1 — 2 /3 Pe 1//2 |’ U l^ 1,/3 ( r he)’ 

and we conclude the proof. □ 

In what follows we denote by Cp a generic constant with depends on (3 but is 
independent of p. 

We finish this Section recalling the following estimates that will be useful later on. 

Theorem 3.3. Let K be a parallelogram in M 2 , 7 an edge of K and —1 < (3 < 0. 
There exist a unique lineal and continuous function T : H l ^{K) -+ such that 

if u E C°°{K ) then T(u ) = u and 

ll^(' u )llrr 0 ’' 9 (7) < Cp\\u\\ H i,P(K), 

Proof. We show the case K = Q = (—1, l) 2 . Let u E C°°(Q), by (J3]) we have that 

u{x,y)= y Ci d J?(x)J?(y). 

i>0,j>0 
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We will bound only ||n(x, —l)||#o,/ 3 m with / = (—1,1), for the rest of the edges we can 
proceed similarly. 


u(x,-l) = 5^ 1) = = I+11. 

i>0j>0 *>0 i>0 j> 1 


then, using dHJ) we have that Jg (— 1 ) = 1 and by (J3J) 

2 ho, P( i) = Y MV < (7o) _1 Y ^ (7o) _1 l u l 2 h°.P( Q ) 

i>0,j>0 


i> 0 


and 


\mi™ m = E (E'M^t-i)) 7 < (e M\jf( 

i >0 j >1 


li 


i> 0 j> 1 

< E (E i^iVi 2 ) (E (-Di’trfr'rV 

i>0 j>1 j>1 

using Lemma 12.31 (12.4)1 and (flUl) there exist a positive constant Cg such that 


¥A\ 2 h 


<c/3 

E(E 

i^iV/)(Eo + 1 ) 2 ' s r 1 )7 


*>o i>i 

j>i 

<C- fl 

E(E 

i>0 j>l 

|cdVj 2 )(/ 


E(E 



*>0 j>l 


< Cg 

Y! 

■\ 2 Tjj 2 Ti < Cp\u\ 2 Hl , HQ) . 


*>0,j>l 


Joining / and II we get 


ima;,-l)||ifo,/j(7) < C' /3 ||m|| // i./ 3(( 3 ) . 

since C'°°(Q) is dense in H l '^( K Q) an unique continuous extension there exists and the 
proof concludes. □ 

For u G H 1 ’ i 3 (K) and 7 an edge of K we denote ||w||#o, 0 ( 7 ) = \\T(u)\\ H o,p^y 

Lemma 3.2. Let I = (—1,1) and a,/3 > —1. Then, there exist C\ t g and C 2 , a such that 
for all polynomials P e V P (I) 


P(x) 2 ( 1 — x^dx < Ci t gp 2 / P{x) 2 ( 1 


P'(x) 2 ( 1 - x 2 ) a+l dx < C 2 , a p 2 / P(x) 2 ( 1 


— x 2 ) /3+1 dx , 

- x 2 ) Q Pr. 
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If, in addition, — 1 < a < ao an d — 1/2 < (3 < —1/4 we can choose C\ and C 2 = 62 ( 0 : 0 ) 
independent of f3 and a. 

Proof. The first inequality can be found in, e.g. (BM971 IBFOOlj . The proof of the 
second inequality is analogous to the proof of Theorem 3.95 of [Sch98| . Following the 
steps in those demonstrations we can see the dependence of (3 and a in the constants. □ 


Lemma 3.3. Let I = (—1,1), —1</3<1 and P G V P (I), there exists v(x,y) which is 
defined on Q with the following properties: 

i) v(x, —1) = P(x)( 1 — x 2 Y, v\dQ\t — 0 where i = I x { — 1}; 

ii) ll^ll// 0 '-/ 3 ^) < Cp(p + l) /3_1 ||P|| // °^(7); 
m) II'^IIh i '-/ 3 (q) < Cp(p + i) /3 ||p|| h °,/3( 7 ). 

If, in addition, 1/2 < /3 < 3/4 we can choose C independent of (3. 

Proof. By Le mm a l3Tl we know that there exists g G V P (I) such that g( 1) = 0, g(— 1) = 
1 and \\g\\ H o,-f 3 ^ < Cg(p + We define v(x,y) = P(x)(l — x 2 )Pg(y), these 

extension obviously satisfies condition i). To prove condition ii) we observe that 


[ v(x,y) 2 (l- x 2 ) ^(l -y 2 ) 0 = [P(x) 2 (l -x 2 fdx [g(y) 2 (l-y 2 ) 0 dy 

JQ Jl Jl 

= miWillsIlF-om < c^(p + i)- 2 < 1 -«||f>|| 2 i ,„,» (/) , 

Hence, 

\\v\\H°-P(Q)<Cp(p + l) (1 /3 ^||P||^o,/3(/). 

To prove condition iii) we calculate ||ff|| h°-P(jQ) and |||4|| # 0 ,-/ 3 ^) 

F)i) / f) \ 

QZ&y) = ^-(x)(l-x 2 y + P(x)^(l-x 2 ) l 3 - 1 (-2x)jg(y), 

then 

J Q {^r(x,y)) 2 (1 - x 2 ) l ^ p (1 - y 2 Y P < C{ (^(x)) 2 (l - x 2 f +1 (l - y 2 )~dg(y) 

+ [ P(x) 2 (l -x 2 Y~ l (l -y 2 )~Pg(y)}, 

JQ 

by inverse estimates of Lemma 13.21 we have that 


(?j?(£)) 2 (1 -x 2 Y +1 dx < Cp(p+ l) 2 J^P(x) 2 (l 
[ P(x) 2 (l — d^Y^dx < Cg(p + l) 2 [ P(x) 2 ( 1 


— x 2 ydx 

— x 2 Ydx. 
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Therefore 


/ (|?(^J')) 2 ( 1 - £2 Y \ l ~y 2 ) 0 < Cp(p+ 1) 2 \\P\\hW(i) 

J Q 


2 


< Cp{p + l) 2 (p + 1) 2(1 i3) \\P\\h°,P(i) 

— Cp(p + 1)' /3 ||P||^o,/3( 7 ). 


On the other hand, 

= P(x)(l - x 2 )^(y), 

then 

J {^(x,y)) 2 (l-y 2 ) 1 ~ l 3 (l-x 2 )-^ = J^P(x) 2 (l-x 2 ydx jT (! -y 2 Y~ P dy, 

by inverse estimates of Lemma 13.21 

^|l(y) 2 ( 1 -y 2 ) 1- ^ < Cp(p + 1) 2 J^g(y) 2 (l-y 2 )~ p dy, 

thus 

< C*/3(P + 1) 2/3 ||-P|llf0,/?(/)• 

from which we conclude the proof. □ 


4 A posteriori error estimation 

In this section we introduce an a posteriori error indicator of the residual type for 
the classical Poisson model problem and, by using the p-interpolation error estimates 
obtained in the previous sections, we show the equivalence between the indicator and 
the error in an appropriate Jacobi-weighted norm up to higher order terms. 

4.1 Problem Statement 

Let 12 C R 2 an open polygonal domain, T = <912 and / € L 2 (12). We consider the 
classical Poisson problem: Find a function u such that: 

f —A u = f in 12 

{ n r (V) 

u = 0 on f 
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The Variational Problem associated to (fTT|) is: Find u G Hq(Q) such that: 

a(u,v) = L(v) V v G Hq(Q) (18) 

where a(u, v ) = f Q Vu ■ Vv and L(v) = f n fv. 

Let T be an admissible partition of D in parallelogram elements, and Sq (T) as in 
m, the Discrete Variational Problem is defined as follows: Find u G Sq (T) such that: 

a(u,v) — L(v ) V v G Sq (T). (19) 

Following the ideas introduced in [ Guo05 |. let Q = (—1, l ) 2 the reference domain, 
(3 > —1, a = (ai,a 2 ) with a* > 0 integer, we define the weight function W^ a in Q as 
follows: 

Wp ta {x i,x 2 ) = (1 -x\Y~ ai {l -xlf~ a2 . 

Note that Wp !<x = Wp- a . 

Then, the weighted Sobolev space H k 'P(Q) is defined as the closure of the C°°(Q) 
function with the norm 

INIl^(Q) = J2 \ \&*v\ 2 Wf), a , 

\a\<k jQ 

by \ u \\k,p(Q) we denote the semi-norm 

MW(Q) = Z) / 1^1 2 ^.«- 

| a\=k JQ 


Let K be a parallelogram and let F : Q —)■ K be an affine transformation. Given 
u G C°°(K) we can define u = u o F G C°°(Q) and (see, for example [Guo05] ) 


\\ U \\H k ’P(K) — \\ u \\H k ’P(Q)- 

Then, for T an admissible partition of fl in parallelogram elements and u G G7 00 (V2), 
we define 


\u\ 


H k ’P(T) 


Ei 

KeT 


\u\ 


H k ’f>(Ky 


Hence, the Jacobi-weighted spaces H k ’@(T) and Hq(T) for T a partition of f1 are 
defined as follows: 

H k 'P(T) =C^0' Kk ’ Pm 

H k /(T) = C' 0 oo (Q) IMI " M(r) . 

Lemma 4.1. Let 0 < /3 < 1 and u G H 1+S (Q) with s > then, u G iL 1,/3 (T). 
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Proof. Since the functions C' 0 O (h2) are dense in H 1 +S (Q), the result can be obtained 
from the inequality (1.8) of [BBOlj on Q, i.e., 

V/effA(Q) ° <A<1/2 (20) 

where p((x,y),dQ ) denotes the distance from (x,y) to the boundary of Q. 

Indeed, there exist ip n G C°°(f2) such that p n —> u as n —> oo in H 1 +S (Q). 

Let K G T and Fk : Q —> K be an affine transformation, we note ip n — u = (ip n —u)oFK- 
Since /3 > 0 


{<pZ~ r u) 2 { 1 - f 2 ) /3 (l - I / 2 )' 3 < / —* 0 as n ->■ oo, 


'<2 


'Q 


then 


||^n — m|| -t 0 as 77 . —^ oo. 

Since 0 < (3 < 1 and (1 — x 2 ) > p((x, y), dQ) we get 


>Q 


(di{<£^u)) 2 {l -x 2 Y 1 {l-y 2 Y< 


Iq p((x,y),dQ) 2X 


with A = yp. Now, since 0<A<l/2by (120|) we follow that 


(dx(<pY-u)Y 


— II 9 x (p n u )\\h x (Q) — \P n u \H 1 + x 


(Q) 


Jq p((x,y), 8 Q) 2X 
then, if A < s we have 

\<Pn - u\ H 1 +A(Q) —> o as n oo, 

and therefore 

[ (9 £ ((^T-«)) 2 ( 1 - x 2 ) /3_1 (l - y 2 Y —» 0 as n G oo. 

JQ 

We can proceed analogously for the other derivative and conclude that ip n —> u as 
n —> oo in H 1,l3 (K). Thus, the result follows since we can do this for any K G T■ □ 

Let u be the solution of (jT8|) and let be the solution of (fT9lh As in [Guo05j . we 
introduce the norm for the error e = u — un denoted by |||e||| and 11 |e| 11 /<- by 


and 


|e||| A '= sup |a(e,u)| < ||e||jj lli8(A:) 
Ihllfii ,-P(K) =1 


|e||| = sup |a(e,n)| < ||e||^-i, 


IHI H 0 w cn =1 


»(TV 
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Let 0 < /3 < 1 and u the solution of (fl 8 ]h if u G H 1+S (£l) with s > AJ1 then 
e G H 1,/3 (T) and, |||e|||# and |||e||| are well defined. 

By (fT 8 l) and (fl9|i we can infer that 

a(e,v N ) = 0 \/v N G S£(T). (21) 


For each / G £° we choose a normal vector ni and denote K in and K out the elements 
that share the edge /. We define 


du at] 
dn 


V( u N \ Kout ) • ri£ - V(u N \ Kin ) ■ n e , 


which corresponds to the jump of the normal derivative of un across the edge t. 

In what follows we assume that 1/2 < j3 < 1 and u, the solution of (ITSlh is such that 
u G H 1+s ( fl) D Hq(Q.) with s > Then, for v G Hq~ i3 (T) such that IMI# 1 - 0 ( 7 -) = 1 

we have that e G i/ 1 ,/ 3 (T) and 


a(e,v) = Y 
K&T 


[ (/ + A%)r + J Y [ 

•* K Z ecdKn£° 


8 un 

dn 



( 22 ) 


For e > 0 and v e G such that ||u — v e \\ H i — e ’ by Theorem 13.31 ||u — 

v e \\ H o-P(e) < Cpe for all £ G £. Since —1 < —j3 < —1/2, if we take Iv e be as in Theorem 
13.21 by using f[ 2 Tj) and (122] ) we get 


a(e, v) = a(e, v — Try) + a(e, Try) 

Y ([ (f+ AuN )( v ~ Iv t) + l E [ 

KgT J K z t>rF>KnF° J£ 


e<zdKn£° 


du 


N 


dn 


(v - Tty)) 


Y ( / (f + Au tv) ((v - V e ) + (v e - Iv e )) 

((v - v e ) + (v e - Tty))) 


Ker 

+ 5 £ /' 

e<zdKn£° 

— E (ll/+ Au tv||HO.0(i<:)(||'tt — u e ||#o,-/3(#) + ||u e —/u e ||#o,-/3(#)) 


KeT 


b E 



dujy 



dn 

1 




V — V e 


I + 11^ — T?y||#o,-/3(^)) • 


25 
































Hence, by (fl5|l . (TT6|) . we obtain 


_ , Q 

a(e,n) < ^ [\\f + Au N \\ H o,p {k) (e + + 2/j Pg (3/2 ~ /3) ||^||g 1 ^(r|^)) 


KeT 

1 
2 


E 


fcafcnf 0 



du N 



dn 

i 


( C P e + _1 , 2 B Pl 1/2 |l' ye ll^ 1 '' /3 ( r l^))) 


ho,P ie)'- ^ -1 + 2/3 

C 


- (ll/ + ^' u Ar||i?o^( jf!r )(e+ + 2 /3 Pk(3/2 ^^\ v ~ v ^\h 1 ,-^t\ uk ) + IMIh 1 --+7+ x ))) 


^er 



duN 



dn 

e 


4 E 

l<zdKn£° 

— (11/ + ^ u n\\h°,P(k) ( e + 


C \ 

[ C P e + _1 I 2/5^ 1/ ^H W ~~ V *\\H',-PiT\ U( ) + W v \\Hi,-P(T\u, e )))j 


ho,p^)' " -1 + 2/3 


c „ (3/2 /3) __ , <E m -(3/2-^) ILill ^ 

)Pk 6 ^-, , o«^A- IMIfP-^TU+J 


Ker 


5 E 


-1 + 2/3 


-1 + 2/3 



du N 



dn 

i 


C' 


(C'+ + Cm, 7 e + 

^o,/3 W V p wt -1 + 2/3 


Pi 1/2 II u IIh i .-/ 3 (ti^)))- 


Since this holds for all e > 0 we conclude that 


a(e,n) < ' ^ (\\f+Au N \\ H o,p( K) p K {3/2 ^ IMI^i, -P{t\ uk )+^ 

+ P KeT V Z ICdKni 

Therefore, 

IH e IH - _l i 2/3 ()\f + Au N\\HO,PiK)PK 3/2 0) + 2 



du N 



dn 

e 


ho,p^) 


- 1/2 

Pi 


KeT 


ecdK n£° 


du 


N 


dn 


ho,p^) 


Pi 


- 1/2 


Let f p = nj , K f as in Corollary 12.21 then 


IPK Pk,K 

C 

~+ 2/3 


e lll < ZY+TTg (ll fpK + AuNWHO^iKiPK^^ ^ + II f ~ fp K \\HO,PiK)PK 


-(3/2 -13) 


KeT 


E 


du 


N 


dn 


H°,p^) 


Pi 


- 1/2 


lcdKn£° 

Let 0 < S < 1/2, we take /3 = 1/2 + 6. Then 
C 

KeT 


- ( WfpK + AunWhO-PWPk 1 5) + 11/ “ fpK \\h°,P(K)Pk 

° KeT 

1 ^ dujy 

9 i fir) HO-P ie.)^' 


-(1-5) 


dn 0 

e<ZdKr\£° U- * 


- 1/2 
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C independent of 5. 

For each element K G T the local error indicator is defined as: 


2 2,2 
Vk = Vb k + Ve k i 

where 

Vb k = 11 f PK + Au n\\ho,p {k) Pk 2 and 

with 

Ve — Il^llif 0 ,/ 3 (;)P^ , Re. = 

Therefore, the global estimator is given by 

V = 2^ Vk- 

KeT 

Thus, we have the following theorem which proves the reliability of the error esti¬ 
mator up to higher order terms. 

Theorem 4.1. Let /3 = 1/2 + 5 with 0 < 8 < \, u the solution of ITR) . un the solution 
of KTO 1) and e = u — un■ Let 77 be as in K dbj) . Assume that u G H 1+S (Q) n Hq(Q) with 
s > ■ Then, there exists a positive constant C such that 

|||e||| < jmax{pf nax ,l}|r/+ ( ^ Pf?\\f ~ fp K \\ 2 H o^ {K )) ' }, 

KeT 

where p m ax = max{ 77 ^| K G T}. The constant C is independent of p and S. 

In order to guarantee that the error indicator is efficient to guide an adaptive re¬ 
finement scheme, our next goal is to prove that px is bounded by the weighted norm 
of the error up to higher order terms. 

The following lemma provides an upper estimate for the first term in the definition 

of (|23D- 

Theorem 4.2. Let /3 = 1/2 + 5 with 0 < 5 < u the solution of i f TR) . un the solution 
of (I22|i and e = u — un■ Let ps K as in Assume that u G H 1+S (Q) D with 

s > Then there exists a constant C(K) such that 

Vb k < C(K) (|||e|||x H-1| f PK ~ f\\HO.P(K))- 

Pk 

The constant C(K ) is not dependent on px and 5. 



(23) 

2 1 \ ^ 2 
ve k = 4 2 ^ vi, 

(24) 

ecdKn£° 


3un 

dn 0 ' 

(25) 


(26) 
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Proof. 


WfpK + ^ u n\\h°’P(K) — II fpK + ^ u n\\ 2 ho,P{Q) 


= / (fpK + ^u N ) 2 (l-x 2 Yil-y 2 Y 

JQ 

= / ( fp K + Au n )vk 

JQ 

where v'k = (f PK + Aujv)(1 — x 2 ) /3 (l — y 2 ) /3 , then we define 

vk = v k ° F k 1 in K and vk = 0 in Q \ K. 

Therefore 

11 J~Pk + ^ u n\\ 2 ho,P(k) = / ( fp K + Aun)vk 

J K 

fv K + / A u N v K + / (f PK - _/>*-). 

JA' J I< 

Now, since vk £ -^o(^) and u solution of (ITS]) then J K fvx = j K V« • Wvk and 



K 


II f PK + Au n \\ 2 h o,p( K) = C(K ) ( J X7u ■ X7v k - J X7u N ■ S7v K + j (f PK - f)v K ^j 
= C(K) ( Ve • Vux + J^(f PK - /K) 

= C'(iT)(a(e,'y if ) + f (f PK - 


= C , (A-')(||nx||Hr-/3(ic)a(e, 




’) + / (fra ~ 


\\ v K\\H 1 -f i (K) JK 

< C(K)(\\v K \\ H i,-f}( K f\\e\\\ K + ||/ PK — /|| J H-o,/3( X )||t>ic||^o,-/3( JS ')). 


Observe that 


I 112 || * 112 

\ v k\\ H o-P{K) = iT-fdl/fO ,-/3(Q) 


= / ^( 1 -^( 1 -^ 


'Q 


= / (/pK + A “iv) 2 (l-£ 2 ) / 3 (l-J / 2 ) /3 = ll/pif + Aujvll^^). 
JQ 


(27) 


On the other hand, 





We note that 


dvx _ fd{f PK + A u N ) 


dx V dx 
and therefore 


(1 - x 2 Y + (. f PK + Au n )P( 1 - x 2 f 1 (-2x )) (1 - y 2 f, (28) 


[ UN) ) 2 {i-x 2 Y +1 {i~y 2 Y 

Iq ox \ Jq dx 


+ / Upk + Au n) 2 {l-x 2 ) p 1 {l-y 2 f) 

Jq 

— C(I + II). 


Now, by the estimates given in Lemma [3.21 it follows that 


I < C Pk / ( f P K + Au n ) 2 (1 - x 2 y\ 1 - y 2 ] 

Jq 


— Cp K WfpK “I" Aujy \\‘ho,/ 3 (k) 


and 


U < Cy\ I (f PK + An 7 v) 2 (l - x 2 f{ 1 - y 2 f = Cp 2 K \\f PK + Au N \\ 2 H0 ,p (K) . 


Therefore, 


’Q 


,dv 


<Q 


Analogously 


(^) 2 ( 1 -^ 2 ) 1 -v 2 ) ^ <Cp 2 K \\f PK + Au n \\ 2 h0 , P(k) . 


j (^) 2 (! - y 2 ) 1 ! \ l - x 2 ) 13 < Cp 2 K\\fp K + A u N \\ 2 H o,p {K) . 


Then, 

\\vK\\m-f>(K) < Cp K \\f PK + Au N \\ H o,P( K y 
Hence, we can conclude that 

II fpK + ^Un\\h°’P(K) < C(K ) (pa'IIM II K + II fp K — /II 5 

and the result follows at once. 

Next, we prove an upper estimate for 7%. 


(29) 


□ 
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Theorem 4.3. Let j3 = 1/2 + 5 with 0 < S < u the solution of (QSjj, un the solution 
of ( E| ) and e = u — un. Let £ G £°, K\ and K 2 G T such that £ = K\ n K 2 and rji as in 
LM) . Assume that u G H 1+S (fl) D with s > Then, there exists a constant 

C{£) such that 

Vi < C{£) | |e| | \k-i_p\ + || f PKl ~ fWm^iK^PK 1 * 6 + |||e||U 2 P* 

+ II fpK 2 - f\\H°^(K 2 )PK 2 +5 + IIMIIn ^ Pf )- 


The constant C(£) is independent of p and 5. 

Proof. Let c<+ = K\ U K 2 and v as in Lemma HOI with P(x) = Re(x, —1), p = p£. Let 
Fk i : Q —> K 1 and Fk 2 : Q —y K 2 affine transformations such that F Ki {I x {-!})= L 
We dehne V£ such that +\#•. = v o F K ' , then 

i) (v e \ e )(x, -1) = Rt(x){l - x 2 Y , vi.\ dLOl \e. = 0; 

ii) IMIff°.-/>(T| u/ ) — C(pe + 1)~ (1 ~®\\Re\\ H o,f>(ey, 

iii) IMIr+-+r|^) < C(p£ + l) /3 ||f?£||^°,/3(^). 

Therefore 


\RA\ho,p w = J Re(x) 2 ( 1 -x 2 Ydx = J R(V t = C(£) J^R e v e 
= C{£) ( - / Au N v t - / Vuat • 


= C'(£)( — / AujyVe — / Vuat ■ Vn£ — / AunV£ — / Vu+r • 
' iiCi J K! J K 2 J K 2 

= C(£) ( - f (f PK + Au N )v£ - [ X7 u n -Vv e + [ (f p - f)vt 
v K\ / 






~ {f P K 2 + Au N)ve- V%-V^+ / ( f P K 2 -f) v £+ f v n- 

J k 2 J k 2 j k 2 Ju> t ' 

It is easy to see that vg G //q(c+) then, since u is a solution of (fT8]l we have that 
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L e f v t = L t Vu ■ Vv * and 

\\Re\\ 2 H o,pm = c(£) (- (f p + A u N )vi - / Vu N -Vv e + / (f p - f)v 

V .!K\ w- ' K - 


'K i 


<K i 


Jk 2 

= c(i) 


( fp K , + Au n )v(. - / Vu N ■ Wve + / (fp - f)v e + Vu- Vv e ) 
Jk 2 Jk 2 Ju e ' 


( fp Kl + &u N )v e + I (/ PJCi - f)v e 


'Ki 




PK 1 


(. f P K 2 + &U N )v e 


K 2 


+ / (fpx 2 ~ f) v e + / Ve-VzA 

J K2 UJg 


(30) 


consequently 

\\Rt\\ 2 H° ’(>(£) — c(t)(\\f rKi + ^ u n\\h°’P(Ki) + II fp Kl ~ f\\H°’t i (K 1 )\\ v e.\\HO--P(K 1 ) 

+ WfpK 2 + ^ u n\\hO’P{K 2 )\\ v i\\ho.-P(K 2 ) + || f P K 2 ~ f\\H°’P(K 2 )\\Vi\\H°-P(K 2 ) 

+ IMkw(ru) a ( e > 


\\ve\\w-P(T\u 


From ii) and iii) we get 


ll-RfllA C(£) (\\f PKl + AuN\\ H o,t3( Kl )(pi + 1) 1/2+5 + II fp Kl ~ + 1) 1/2+5 

+ II fpx 2 + ^ u n\\h°’P(k 2 )(Pi + 1) 1/2+5 + \\f PK2 - f\\ H o,P( K2 )(pe + 1) 1/2+5 

+ (p«+ l) 1 / 2 +< ll|e||ln.,). 

then 

Pi / ll^ll A C(£)^\\fp Ki + AunWho^^k^Pz 1+5 + II fp Kl — 1+5 

+ II/pk 2 + ^ u N\\H0-P(K 2 )Pi 1+5 + ||/pk 2 _ f\\HO,P(K 2 )Pl 1+5 + lll e ll|r|^P 5 )- 

Since the polynomial degrees of neighboring elements are comparable, from Theorem 
14.21 we can deduce that 

p 7 1/2 H^II^W < ^(lIMIkiP 5 + II f PKl ~ f\\HW{Ki)P~K\ +6 + lll e llk 2 p 5 

+ ll/pjf 2 _ f\\H°-P(K 2 )PK 2 +S + lll e lllr|^P 5 )) 

and we conclude the proof. □ 

Now we may conclude the efficiency of the error indicator. 
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Theorem 4.4. Let /3 = 1/2 + 5 with 0 < 8 < 1/4, u the solution of [TED , un the 
solution of [T9 j ) . e = u — un and rj as in [2BD . Assume that u e H 1+S (Q) D with 

s > . Then there exists a constant C independent of p and 5 such that 

V < Cmax{p s m a X ,l}{\\\e\\\ + (^PlfWfpK ~ fllm^R)) 

KeT 

Proof. Let K e T, and Vk as in the proof of Theorem 14.21 in this proof we show that 
II f PK + Awtv \\ 2 ho.P( K) = C(K)a(e, v K ) + C{K) j (f PK - f)v K , 

J j 

then, 

Vb k = a (+ (Pk + 1)~ 2 C(K)v k ) + (p K + l)- 2 C(K) [ (f PK - f)v K , 

Jk 

let Vj- = J2kgt(Pk + l) -2 C , (-^0+fo from (1271) and (l29|i we find that 


Vbk = a(e, V T ) +^(PK + 1 ) 2 C{K) j (f PK - f)v K 

KeT KeT 


hence, 


< 


v T 


ll^rll/p-+r) 
e||| ^2 ( P K + ^)- 2 C(K)\\v K \\m-P(K) 


) + X>* + / (f PK -f)vK 

K&T ^ K 


KeT 


+ ^(PX + 1) 2 C(K)\\f pK — f\\H 0 ’P(K)\\ v K\\H°-P(K) 
KgT 

A C , (lll e lll ^ 1 |l fpK + Au N \\ H o,P(K) 

Ker 

+ ^2(pk + 1) _2 ||/ Pk - f\\HW{K)\\fp K + Anjvll H°’P(K)) 


KeT 


< 


C (lll e lll + ( ^ ( Pk + 2 \\fpK - f\\HO'f>(K)) / ) ( Y. (pk + 1) 2 ||/pjf + ^ u n\\ho,P{K)) 


1/2 


A'eT 


KeT 


Y Prk) / A ^(jIMII + ( Yj^ )k + 1) 2 \\fp K + ^ u n\\ 2 ho,p(k)) ^ J • (31) 

KeT KeT 

Let £ G £°, and V£ as in the proof of Theorem 14.31 from (1301) we have that 

\\Rt\\ 2 H w{i) = °( e > C{£)vi) + C(£) ( - f ( f PKi + A u N )vt + j ( f PKi - f)v e 


Ki 


IK i 


(f P K 2 + A u N )v e + / {f PKo - 


'k 2 


'k 2 


PK 2 


32 






Then, 


*52 Vt = ||^C'(£)p/^||^w( r) a(e, 


ee£ 


eee 


Etef cWp<V 


+ ^ ( - / (fpK 1 + Aujv)^- 1 + [ (f PK - f^ePi 1 

ee£ Jk ' 7a 'i 

- / (/wr 2 + A^jv )^ -1 + / Upk 2 ~ f) v ePl) 

-I K2 J K'r. ' 


>k 2 


K 


l e lll y ^C(£)p e 1 \\v£\\ H i,-p(r\ Ui ) + ll/pjg, + AmjvIIho,^^) 


fe£ 




+ \\fp Kl ~ fWHO’^K!) + \\f P K 2 + Au N \\ H o,0(K 2 ) + \\fp K2 - f\\H 0 'P(K2))\\ v e\\H 0 ’-P(T\u e )Pl 1 - 


i LU £' 


From the estimates for vg, given in ii) and iii) in the proof of Theorem 14.31 it follows 
that 


J2^ A C\\\e\\\J2Pe l 2 Pe\\ R i\\HO^(i) + J2 v\fp Kl 


ee£ 


ie£ 


A Ur 


£G£ 


N\\H°>P(K\) 

+ \\fpK 1 - f\\H°’f>(Ki) + II fpK 2 + ^Un\\h°’P(K 2 ) + II f P K 2 - f\\H°’P(K 2 )) II^IU°^(^)Pf 1/2 Pe 1+5 

< CpLJMlCErf) 1 ' 2 + c£p( ( ‘ 1 +<) (ll./W, + A« N |li..», Kl) + II/«, - ff m , (Kl) 


£&£ 


£e£ 


+11 fpx 2 +Au N \\ 2 H o,p( K2) +11 f PK2 ~ f\\ 2 H °.p( K2 )) / iy^.vi) f ■ 


2 \l /2 


ee£ 


Since the polynomial degrees of neighboring elements are comparable we have that 

\ 1 / 2 ' 


QZ^ 2 ) 1/2 - C {Pm<Ell e lll + {^2pkVk+Pk L+d) \\fpK - /Ilk/3(X)) 


( 5Z ^) V2 + (J^PkUpk - f\\ 2 HW{K)) 


KeT 


£&£ K&T 

Then, 

(£.I ) 1/2 < cp s max { 

£e£ 

and by (j3T)) 

(5Z^ 2 ) V2 A Cp. 

£&£ ' K&T 

Hence, we are in conditions to compute 77 


1/2 


k&t 


/ 

max 


£ Pi^ll/E f\\ 2 H°’P(K)) 


1/2 


V 2 = ^2(.Vb k +Ve k ) < C ( £ +£ 7 7<) - C ' max {^a^ 1 }{lll e lll + ( ^PxWfpK -/ll^(iC)) 1 


KeT 
as claimed. 


KeT 


ee£ 


KeT 


□ 
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We end the paper by emphasizing that, as far as we know, the quasi-optimal esti¬ 
mates reached in Theorems 14.11 and 14.41 are the best results that can be obtained for 
error estimators of the residual type for the two dimensional case. 
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